Vacuum Energy and Repulsive Casimir Forces in Quantum Star Graphs 
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Casimir pistons are models in which finite Casimir forces can be calculated without any suspect 
renormalizations. It has been suggested that such forces are always attractive, but we present 
several counterexamples, notably a simple type of quantum graph in which the sign of the force 
depends upon the number of edges. We also show that Casimir forces in quantum graphs can be 
reliably computed by summing over the classical orbits, and study the rate of convergence of the 
periodic orbit expansion. In generic situations where no analytic expression is available, the sign 
and approximate magnitude of Casimir forces can often be obtained using only the shortest classical 
orbits. 

PACS numbers: 03.70.+k, ll.10.Kk, 42.25.Gy, 03.65.Sq 



I. INTRODUCTION 

According to a classic calculation [l[ , the Casimir force 
inside a roughly cubical rectangular parallelepiped is re- 
pulsive; that is, it tends to expand the box. The rea- 
soning leading to this conclusion is open to criticism on 
two related grounds: It ignores the possibility of non- 
trivial vacuum energy in the region outside the box, and 
it involves "renormalization" in the sense of discarding 
divergent terms associated with the boundary although 
(unlike the case of parallel plates, or any calculation of 
forces between rigid bodies) the geometry of the bound- 
ary depends upon the dimensions of the box. Recently 
(see also 0) a class of scenarios called "Casimir pistons" 
has been introduced to which these objections do not ap- 
ply. The piston is an idealized plate that is free to move 
along a rectangular shaft, whose length, L~a, to the right 
of the piston is taken arbitrarily large (Fig.[T]). Both the 
external region and the divergent (or cutoff-dependent) 
terms in the internal vacuum energy are independent of 
the piston position, a, so that a well-defined, finite force 
on the piston is calculated. One finds that this force is 
always attractive, both for a two-dimensional scalar-field 




FIG. 1: A rectangular piston in two dimensions (cf. 3). In 
three dimensions there is another length, 62 , perpendicular 
to the plane of the figure. 
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model with the Dirichlet boundary condition Q and for a 
three-dimensional electromagnetic field with the perfect- 
conductor boundary condition [J]. 

Barton Q showed that the piston force can be repul- 
sive for some (not too small) values of a if the conduct- 
ing material is replaced by a weakly polarizable dielec- 
tric. This result is somewhat ironic in that one reason 
for suspicion of repulsive Casimir forces is the belief that 
the force between disjoint bodies of realistically modeled 
material should be always attractive. The unexpected 
result is easily understood, however, as being due to at- 
traction between the piston and the distant part of the 
shaft. The effect would disappear if the shaft extended 
a long distance to the left of the fixed plate ( "baffle" ) at 
a = as well as to the right of the piston. 

In the present paper we study the vacuum energy and 
Casimir forces in one-dimensional quantum graph mod- 
els and observe several situations with idealized bound- 
ary conditions for which the piston force is unambigu- 
ously repulsive. In quantum graphs of high symmetry, 
the Casimir forces may be calculated analytically. More 
generally, we show that these forces may be obtained sys- 
tematically from a sum over the classical periodic orbits 
in the graph, as done in three-dimensional problems in 
[H 0, Hf, and we discuss the rate of convergence of the 
periodic-orbit expansion. In some cases, the sign and 
approximate magnitude of the force on a Casimir pis- 
ton may be obtained using only the shortest orbit hit- 
ting that piston. Although the quantum graph models 
are less realistic than those studied in [| and @, they 
do show that repulsive Casimir forces do arise physically 
and are not inevitably an artifact of a naive renormaliza- 
tion scheme. Our effects are unrelated to that in [5( and 
do not depend on the asymmetry noted above in con- 
nection with that paper. The periodic-orbit techniques 
discussed here have relevance to the study of Casimir en- 
ergies in more realistic geometries (cf. §), including two- 
and three-dimensional chaotic billiards. In an Appendix, 
we consider a situation in which an unambiguously repul- 
sive Casimir force appears for the electromagnetic field 
in a three-dimensional geometry 
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Throughout, we take Ti = 1 



c. 



II. VACUUM ENERGY IN QUANTUM 
GRAPHS 

A finite quantum graph [l(| El, E3, EH consists of 
5 one-dimensional undirected bonds or edges of length 
Lj (j = 1,...,B). Either end of each bond ends at 
one of V vertices, and the valence v a > 1 of a ver- 
tex is defined as the number of bonds meeting there. 
A normal mode u of the quantum graph has the form 
Uj(x) — o,jcos(kxj) + bjSm(kxj) on every bond j, and 
satisfies the specified boundary conditions at each ver- 
tex. Despite their simplicity, quantum graph models have 
previously shed light on a number of important physical 
problems, having served originally as models of conju- 
gated molecules, and more recently of quantum, electro- 
magnetic, and acoustic waveguides and circuits. These 
models have also served as valuable testing grounds for 
studying more general properties of quantum behavior, 
including Anderson localization, quantum chaos, adia- 
batic quantum transport, and scattering. A recent re- 
view may be found in [14j |. 

In the spirit of abstract modeling, the vacuum en- 
ergy of a graph is defined as the sum (renormalized) of 
zero-point energies over all normal-mode frequencies cu n , 
where the frequency uj n is equal to the wave number k n in 
our units. It is convenient to apply an exponential ultra- 
violet regularization (the same answer would be obtained, 
for example, by a calculation with zeta functions): 
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where 



(1) 



(2) 



n=0 



is the trace of the so-called cylinder kernel [If 



III. ANALYTIC EXAMPLES OF REPULSIVE 
CASIMIR FORCES 



One-dimensional piston with mixed boundary 
conditions 



The first example is already rather well known, in its 
essence. Consider a scalar field quantized on a line di- 
vided into three parts by two points, at each of which 
either a Dirichlet or a Neumann boundary condition is 
imposed. The contributions of the two infinite (or, bet- 
ter, extremely long) intervals to the Casimir force will 
vanish. (As emphasized in [ij , the force contributed by a 
long shaft is entirely associated with periodic orbits per- 
pendicular to the shaft, which do not exist in the one- 
dimensional case.) Let the length of the central interval 



(3) 



be a. Then the frequencies of the normal modes are 

rm 

u n = — 
a 

for nonnegative (or positive) integer n, if the boundaries 
are both Neumann (or both Dirichlet, respectively), and 
one has 

oc 

T(t) = Y e- nnt/a 
n=Q,l 
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Thus the regularized vacuum energy is 



E(t) = 



27rf 2 



2^ + °W 



(5) 



The leading, divergent term is proportional to the inter- 
val length a and corresponds to a geometry-independent 
constant energy density. This term is compensated in 
the force by similar terms in the exterior regions, already 
discarded. Then letting t — > 0, we obtain the well-known 
attractive force 

BE TT 

F ^-lTa=-2^- (6) 
More precisely, if the entire space has length L, then 
the regularized energy of the exterior regions is 

L 



27rf 2 



+ 0(£- 1 ). 



(7) 



The second term is negligible as L — > oo, and the first 
term combines with the first term of ([5]) to make a term 
independent of a, which, therefore, is an unobservable 
constant energy shift that contributes nothing to the 
force. Henceforth we shall not repeat this type of ar- 
gument every time it is needed, and will simply refer to 
such endpoints as Neumann or Dirichlet pistons. 

On the other hand, if one boundary is Dirichlet and 
the other Neumann, then the eigenfrequencies are 



(2n + 1)tt 
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The same calculation leads to 



T(t) = e 
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the regularized energy is 
E(t) = 



27rf 2 



and the force comes out to be repulsive: 

TT 



F = 



48a 2 



(10) 
(11) 
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2a, respectively. When B > 3, however, the force is re- 
pulsive: if the pistons are free to all move together, they 
will tend to move outward. (More generally, a periodic- 
orbit calculation, such as discussed in Section ITVl is ap- 
plicable to unequal piston displacements and confirms 
that the force on each individual piston is outward, so 
there are no other, asymmetrical modes that are partly 
attractive.) This repulsive effect cannot be attributed to 
mixed boundary conditions, since all the conditions are 
of the Neumann type. (However, replacing all the pis- 
tons with Dirichlet pistons while maintaining Q12p would 
interchange the roles of the two types of eigenvalues and 
produce attraction for all i? > 1.) 



FIG. 2: A star graph with a piston installed in each edge. 
(The pistons are actually points; the edges have no thickness.) 



IV. PERIODIC-ORBIT CALCULATIONS FOR 
GENERAL GRAPHS 



B. Quantum star graphs 

In the next model the space consists of B onc- 
dimensional rays of large length L attached to a central 
vertex (Fig. [2]). In each ray a Neumann piston is located a 
distance a from the vertex, so that a normal mode of the 
field in ray j must take the form Uj(x) = Cj cos(ui(x — a)) 
when x is measured from the center. At the central ver- 
tex the field has the Kirchhoff (generalized Neumann) 
behavior 



Uj (Q) = C for all j, 
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(12) 



The following analysis is part of a broader stud y o f vac- 
uum energy in quantum graphs [lj| (see also [I^ll8l[l9l|). 

There are two types of normal modes. First, if 
cos(a-ia) 7^ 0, we have from (fT2"|) that Cj = C / cos(wa) 
and tan(wa) = 0, whence uj is one of the numbers ([3]). 
Second, if cos(wa) = 0, then id is one of the numbers ([8]) 
and 



£9 = 0, 



(13) 



which has B — 1 independent solutions. Therefore, the 
energies and forces are just the appropriate linear combi- 
nations of those calculated in the previous example: the 
regularized energy for the whole system is 



E(t) 



BL (B - 3)tt 



+ 0{L- 1 ) + O(t) , (14) 
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and the force (either from (|14p or from © and (fTTjO is 
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When B = 1 or B = 2, the result reduces properly 
to that for an ordinary Neumann interval of length a or 



For a general quantum graph, e.g., for a star graph 
with unequal bond lengths or with more complicated 
boundary conditions, no simple expressions for the 
normal-mode frequencies ui n are available, and thus the 
vacuum energy and Casimir forces cannot be computed 
in closed form. Computing the spectrum numerically, as 
discussed below, allows for an accurate evaluation of the 
vacuum energy for any specific quantum graph, but this 
type of brute force calculation must be repeated anew 
for every geometry and does not provide much physical 
insight regarding the attractive or repulsive character of 
Casimir forces in different cases. Instead, much intu- 
ition may be obtained using a classical-orbit approach, 
where the sign and magnitude of every contribution to 
the vacuum energy are seen to be directly related to bond 
lengths and boundary conditions at the vertices. 

It is convenient to describe boundary conditions at ev- 
ery vertex a by a unitary v a x v a scattering matrix a a 
(which acts on the space of undirected bonds meeting at 
vertex a). For example, a Neumann or Dirichlet bound- 
ary condition at a vertex of valence v a = 1 corresponds 
to a scattering matrix a a = (+1) or (—1), respectively, 
while the Kirchhoff boundary condition is described by 
(cr a )jj' = — — Sjj>. Together these constitute a 2B x 2B 
scattering matrix S for the entire graph of 2B directed 
bonds [HI EH or bond-ends [H E0, Q . To make the 
following arguments valid, we must assume that S is in- 
dependent of energy or frequency (^-independent), as is 
true for the Dirichlet, Neumann, and Kirchhoff bound- 
ary conditions we treat here (but not for the more general 
Kirchhoff-type boundary conditions where a potential is 
attached to each vertexjSl] and p], E El))- Then 
one can construct [Til . [l3| a trace formula relating the 
spectrum of a graph (away from the point u> = 0, which 
makes no contribution to vacuum energy anyway) to its 
periodic orbits, 
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(16) 



p r—1 



(Variations on the trace formula have been found in [Iq1 |, 
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dUi 0! an d elsewhere.) In (fl~6f the values uj n are 
the normal-mode frequencies, L = Ylj=i Lj is the to- 
tal length of the graph, which determines the smooth 
(Weyl) contribution to the spectrum, and the sum over p 
is a sum over primitive periodic orbits (orbits that can- 
not be written as repetitions of shorter orbits). Each 
p takes the form p = jxj% ■ ■ ■ j n where every ji is a di- 
rected bond. The corresponding amplitude of the prim- 
itive periodic orbit is given by a product of scattering 
factors, A p — Sj 1 j 2 ■ ■ ■ Sj n _ 1 j n Sj n j 1 , the metric length of 
the primitive orbit is L p = Lj 1 + • • • + Lj n , and each r is 
a different repetition number of our base primitive orbit. 

Substituting the spectrum given by Eq. (fT6|) into 
Eq. fTj), we obtain 



E{t) = 



27rf 2 



r 4eet4+ w > 

2tt * — ' ^ — ' Lr>r z 



p r—1 



(17) 



As discussed previously, the finite vacuum energy, which 
is relevant for computation of Casimir forces, is obtained 
by dropping the divergent Weyl term and taking the limit 
t -> 0, 



1 oo 



p r—1 



L n r 2 



(18) 



A mathematically rigorous derivation and proof of (con- 
ditional) convergence of Eq. (fT8)) will appear in [19( • 

Equivalently, we may begin with the free cylinder ker- 
nel in one dimension, 



t 



T (x,x',i) = - 



1 



7T (X - X 1 ) 2 + t 2 



(19) 



apply the method of images to include scattering from 
the vertices, take the trace 



T(t) 



dx T(x, x, t) 



^+ Re EE- 



— TT {tL p ) 2 
p r—1 v y/ 



O(t 2 ),(20) 



and finally use Eq. (JTJ) to obtain the result (fTTj) . (This 
construction, which generalizes the study of the heat ker- 
nel in is described in detail in [HI].) 

The Casimir force on any piston may be obtained easily 
by differentiating Eq. term by term with respect to 
the appropriate bond length Lj . 

We note that the expansion (fT8|) of the vacuum energy 
is exact and involves periodic orbits only. The derivation 
of Eq. ([20]) hinges on the identity (a a ) 2 — I for the scat- 
tering matrix at each vertex. This condition holds for any 
k- independent scattering matrix [20 . [2l| , including real 
scattering matrices of the form used here, but also com- 
plex energy-independent scattering matrices in the case 
of time-reversal symmetry breaking by magnetic fields. It 
is the crucial ingredient in proving that closed but nonpe- 
riodic paths (i.e., paths that start and end at x but with 



opposite momenta) make no net contribution to the vac- 
uum energy. When S depends on k, two complications 
arise. First, the method of images cannot be so easily ap- 
plied to "time-domain" integral kernels such as T and the 
heat kernel, because the reflection law becomes nonlocal 
in t. Second, the identity (c^) 2 = I no longer applies, 
and the nonperiodic paths make a nontrivial contribu- 
tion to the vacuum energy (and to the density of states, 
Eq. (|16p . even when u / 0). Both effects are visible in 
the investigations of the simplest special cases in [l?], [HJ . 

To evaluate the accuracy of the periodic-orbit expan- 
sion in situations where no analytic expression for the 
vacuum energy is available, we may compare with a 
brute-force calculation where the spectrum is evaluated 
numerically. For a general ^-vertex graph, the normal- 
mode frequencies are given by solutions of a characteristic 
equation det h(u>) — 0, where h(uj) is a V x V matrix [ill ]. 
For the special case of a star graph with irrationally re- 
lated bond lengths, we have 



B 

E 

3=1 



tan(wLj +8j) = Q, 



(21) 



where 9j = or ir for a Neumann or Dirichlet piston on 
bond j, respectively. In any case, given a method for 
obtaining a numerical spectrum oj n , we may evaluate 



^finite (t) = ^2 -ui n ( 



— UJrit 



L 
2^rf2 



(22) 



to any desired accuracy by summing over all aj n < w max 
where u> max l/t. Since the divergent term associ- 
ated with the Weyl density of states, or equivalently with 
the free one-dimensional geometry, has already been sub- 
tracted, we only need take the numerical limit t — > to 
obtain the true vacuum energy E c . Expressing the reg- 
ularized vacuum energy as a power series, 



£finite(i) = E c + ait + a 2 t 2 



(23) 



we may apply Richardson extrapolation to approximate 
the vacuum energy to any desired order of accuracy, E c — 
^numerical + 0(t s ), by evaluating E &nite (t) at s distinct 
values of the regularization parameter t. 



V. RATE OF CONVERGENCE OF 
PERIODIC-ORBIT EXPANSION 

We consider a star graph with Kirchhoff boundary con- 
dition for B bonds meeting at the central vertex, and a 
Dirichlet or Neumann piston on each bond at a distance 
Oj from the central vertex (i.e., the pistons may be lo- 
cated at different distances from the center). The leading 
contribution to the vacuum energy is given by the short- 
est primitive orbits, each of which travels back and forth 
along a single bond. Including all repetitions of such or- 
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FIG. 3: The force on a piston in a star graph with B bonds of 
length 1, Kirchhoff boundary condition at the center, and ei- 
ther Neumann or Dirichlet boundary condition at each piston 
is computed using only the shortest periodic orbit (Eq. (|24|l) 
and compared with the exact answer. Positive values indicate 
repulsive forces. 



bits, we obtain 



e: 



shortest 



1_A™ 1/2 



r cos(r9j 



(24) 



where 8j = for a Neumann piston or ir for a Dirichlet 
piston. For example, for all Neumann pistons the sum 
over r can be evaluated as a dilogarithm, which in turn 
can be expanded in powers of 1 / B as 



E: 



shortest 



vr / 24 In 2 
~ 48 V 7r 2 B 



(25) 



This approximation compares well to the analytic result 
2g (l - |) f for B equal-length bonds (Eq. USD ). 

The results are illustrated in Fig. [3J where the exact 
force on each piston in a star graph having either all 
Dirichlet or all Neumann pistons is compared with the 
contribution to the force from the shortest periodic or- 
bit. The repulsive behavior in the Neumann well 
as the attractive behavior in the Dirichlet case, are well 
explained by considering only the shortest periodic or- 
bit, i.e., the bounce between the piston and the central 
vertex. 

To obtain a better approximation, we may systemati- 
cally include contributions from longer orbits. In Fig. 21 
we show the convergence of the sum ([T5]) when all or- 
bits, including primitive orbits and repetitions, of total 
length rL p < L max are included in the summation. In 
this example, a star graph with B = 4 bonds, all Neu- 
mann pistons, and unequal bond lengths is used, so the 
exact answer is obtained to the necessary accuracy from 
a numerical spectrum as described in Section IIVI We 
note that the rate of convergence is given by 



\Ei 



E r 



1 



(26) 



FIG. 4: The error \E L C 



E c \ in the periodic orbit expan- 



sion for the vacuum energy is shown for a star graph with 
four bonds of length 1.1, 1.6176, 1.2985, and 1.1159, and a 
Neumann piston at the end of each bond. 



consistent with the fact that each contribution to Eq. |T 
from orbits of length rL p G [L max ,L max + A] scales as 
L~^ x for large L max , and all such contributions appear 
preferentially with the same (negative) sign. 

In more general situations, involving non-star topolo- 
gies, more complicated boundary conditions, or non-zero 
gauge fields, orbits of different length are expected to 
contribute with random signs to the sum (|18[) . The er- 
ror made by omitting orbits of length greater than L max 
takes the form Y^=o ® n > where D n , associated with all 
orbits of total length rip £ [L max +nA, L max +(n+l)A], 
scales as D n ~ (L max + nA)~ 2 , but the D n appear with 
random (uncorrelated) signs. The mean squared error 
then scales as J2n=o D l ~ E^Lo^max + nA)" 4 ~ L. 
and the root mean square error decays as 



3 

max ? 



\Ei 



-3/2 



(27) 



As an example, in Fig. [5l we consider the convergence of 
the periodic-orbit sum for the same 4-bond star graph, 
but with a Dirichlet instead of Neumann piston on one 
of the bonds. The behavior is consistent with the faster 
rate of convergence predicted by Eq. (|27p . 



VI. SUMMARY 

We have seen that unambiguously repulsive as well as 
unambiguously attractive Casimir forces arise in simple 
quantum-graph models, and that the sign of the force 
in a given geometry may often be easily understood in 
terms of the short periodic orbits of the system. We 
have also examined (numerically) the rate of convergence 
of the periodic-orbit expansion. Classical-orbit approx- 
imations may also be useful for understanding the sign 
of Casimir forces in higher-dimensional piston systems 
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FIG. 5: The error |£' ( f max —E c \ is shown for the same quantum 
graph as in Fig. |4j but with a Dirichlet piston at the end of 
the first bond. 



where no analytic solution exists, for example, in two- or 
three-dimensional chaotic billiards. 
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APPENDIX A: 



INFINITELY PERMEABLE 
PISTON 



In principle, a repulsive piston can be constructed in 
the more realistic case of the electromagnetic field in di- 
mension 3, in analogy with our original one-dimensional 
model. If the electromagnetic analog of the Dirichlet con- 
dition is a perfect conductor, then the analog of the Neu- 
mann condition is a material with infinite magnetic per- 
meability, and the Casimir force between slabs of these 
two different types is repulsive (25|. (A list of references 
on this topic appears in |2y|.) The existence of real mate- 
rials with sufficient permeability to exhibit Casimir repul- 
sion in the laboratory is controversial [27], [28|, |29( . Here 
we merely check that the piston effect discovered by Cav- 
alcanti Q and the MIT group Q does not destroy the 
repulsion shown by less sophisticated calculations. This 
is not trivial, since the effect arises from the action of the 
shaft walls on the transverse behavior of the field. 

Following Lukosz but in a notation closer to Cav- 
alcanti's (see Fig. [T}, we consider a rectangular box with 
dimensions a, b\ , and 62 ■ As previously exemplified, 



we can calculate a finite vacuum energy naively, in full 
confidence that the discarded divergent terms will can- 
cel when a force is calculated for the piston system as a 
whole. We are interested in the case where the piston 
(the surface that is free to move) is infinitely permeable 
but the shaft and the baffle (the rest of the box) are per- 
fect conductors. By the Rayleigh-Dowker argument [30] ; 
the energy, E a , of such a box is 



E a — E 2a — E a , 



(Al) 



where E a is the energy of a totally conducting box also of 
length a. By differentiation with respect to a (not 2a!), 
this relation extends to forces and pressures. (Through- 
out this discussion "pressure" simply means "force per 
area" without necessarily implying a local pressure inde- 
pendent of position on the wall.) Thus (fTTj) follows from 
© by virtue of 



(A2) 



and the three-dimensional analogs will involve quantities 
proportional to 
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(A3) 



When a -C 6, , Lukosz calculates an attractive pressure 



Pn .= 



which implies by 



240a 4 ' 
1J Boyer's formula 2i 



240a 4 



(A4) 



(A5) 



for the box with one permeable wall. The external (long) 
part of the shaft has length L — a 3> 61 = 62 = b. For 
this limit, Lukosz finds a repulsive pressure (involving 
Catalan's constant) 



P = 



0.915965 
24o 4 



(A6) 



Just as in [H , the resulting force is inversely proportional 
to the cross-sectional area and is independent of L — a, so 
the corresponding energy term is proportional to L — a. 
Therefore, application of (|A1[) gives 



P L -a = P L -a = 



0.915965 
246 4 



(A7) 



(as ought to be the case, since the nature of the plate at 
the distant end of the long shaft ought to be irrelevant). 
To find the total force on the piston, we must reverse the 
sign of (|A7|) . add it to (|A5|) . and multiply by the area, 
b 2 . The point is that the total force is positive if a -C b; 
the long external part of the shaft has negligible effect in 
that case. 
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On the other hand, for a cube Lukosz found that the 
perfectly conducting box was already repulsive. The for- 
mula (|Alj) does not yield a simple factor — \ in that case, 
because the doubled box is no longer a cube. Neverthe- 
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